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On the effect of overshooting as predicted by the modelHng of the 
pre-main sequence evolution of a 2 solar-mass star 
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ABSTRACT 

We discuss the effects of convective overshooting in the PMS evolution of intermediate mass 
stars, by analysing in detail the early evolution towards the main sequence of a 2 AIq stellar 
model. These effects can be extremely important in the end of the PMS, when the abundances 
in CNO elements approach the equilibrium in the centre. We provide a possible physical 
explanation on why a moderate amount of overshooting produces, as the star approaches the 
ZAMS, an extra loop in the evolutionary tracks on the HR diagram. 

An interesting feature is that there is a very well defined amount of overshooting (for a 
given stellar mass and chemical composition) beyond which a loop is produced. For smaller 
amounts of overshooting such a loop does not take place and the evolutionary tracks are 
similar to the ones obtained by Iben ( 1965). The amount of overshooting needed to produce 
the loop decreases with stellar mass. 

We discuss the underlining physical reasons for the behaviour predicted by the evolution 
models and argue that it provides a crucial observational test for convective overshooting in 
the core of intermediate mass stars. 
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1 INTRODUCTION 

Convection and convective overshooting can play a fundamental 
role in the evolution of a star, by being efficient mixing processes. 
From the coupling of convection with nuclear reactions, signif- 
icant changes on the evolution may be expected when convec- 
tive overshoot is present, regardless of using different prescrip- 
tions for convection and/or overshooting. Several authors have ad- 
dressed the effect of including overshoot in mod elling the evolu- 
tion of main and post-main sequence stars (e.g. IStothers & ChinI 
[l99Q; Mowlavi & Forestini 1994, and references therein). But the 
problem of how much is the amount of overshooting that has to 
be included is still far from solved (see ..Renzini ,1987). Several 
authors modelled core overshooting in stars ("e.g. | Roxburghlll978t 
iBressan et al]ll98lHLange3ll98dKion3ll985lll98d ). arriving at 
different results concerning the extent of the convective penetration 
into the radiative zone, from negligible to quite important. Another 
open issue is how oversh oot affects the local temperature stratifi- 
cation te. g. ICanutol 1 997h . But as far as evolution is concerned the 
exten t of overshooting is the key aspect jPeupreellOOCtlRibas et alj 

Eqq3). 



There are several indications, both theoretical and observa- 
tional, that support the existence of overshooting in the convective 
core of intermediate and high mass stars. The comparisons of mass 
and radii of eclipsing binaries with theoretical models suggests not 



only the existence of a certai n amount of overs hooting, but also its 
dependence on stellar mass jRibas et alj|2000l) . On the other hand 
theoretical predictions of the apsidal motion rate for very eccen- 
tric binary orbits can be compatible wit h observed values, only i f 
overshooting is included in the models (Claret & Gimenez 199" 
Moreover 2D hydrodynamic simulations performed by Deuorei 
(2000) for stars with masses from 1.2 to 20 Mq predict convec- 
tive core overshooting for all models. 

The effect of overshoot on stellar evolution has already been 
extensively discussed in the literature, but its effect on models of 
pre-main sequence stars has never been discussed taking the full 
network of nuclear reactions into account. The overall effect to be 
expected is a slight delay in the evolution, redefining the age at 
which the star arrives at the zero age main sequence (ZAMS). As 
the zero point for age of young stars is in itself an open problem, 
such a global effect of overshoot on the PMS has been ignored. 

The mixing at most of the overshoot layer at the border of a 
conv ective core in main se quence stars is expected to be efficient 
(e.g. browning et aljl2004) . This should also be the case for the 
end of the PMS evolution, as the time scales for mixing due to con- 
vection are still much smaller than the time scales for nuclear reac- 
tions. The uncertainties on the thermodynamic stratification within 
an overshoot layer are not expected to be of major relevance for 
determining the effects on the early evolution, but the extent of the 
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overshoot layer and its effect on mixing the stellar material at that 
location are bound to be decisive. 

In a previous work I'Maraues et al."2004') we studied the ef- 
fects of the time step and convective overshooting in the calibra- 
tion of the young binary EK Cep. It has been shown that the time 
step adopted for the PMS evolution is a key numerical aspect to 
adequately reproduce the expected behaviour in this rapid phase of 
evolution. When the time step is properly defined for calculating 
the evolution of intermediate mass stars with overshoot the evo- 
lutionary tracks near the end of the PMS display an extra loop. 
This loop seems to be present in the output of previous works (e.g. 
[siess et al. 2000) although it has never been discussed in the liter- 
ature, as for most cases an inadequate time step in the integration 
procedure would mask the loop as being a numerical artifact. Con- 
sequently its nature and physical validity has never been analysed 
or discussed in the literature. 

Given that an adequate definition of the time step for the PMS 
evolution together with the most recent set of physics produces the 
loop near the ZAMS (when there is a significant amount of over- 
shooting), it has lead us to perform a detailed analysis on the un- 
derlying reasons why such a loop is predicted by an evolutionary 
code. In particular the goal of this work is to establish if such an 
effect should be expected in real stars. Instead of dismissing it as 
a numerical artifact, as implicitly done by other authors, we have 
investigated if it could be consistent with the physical picture that 
the new state-of-the-art physics predicts (in particular the network 
of nuclear reactions). Here we provide such an analysis showing 
that there may be a physical reason for such a behaviour. In spite of 
the difficulties of modelling convection and convective overshoot, 
our analysis seems to indicate that the loop in the HR diagram at 
the end of the PMS evolution could happen in real stars. 

We start by discussing how overshoot is introduced in a one 
dimensional evolutionary code and move on to show evolutionary 
tracks without and with overshoot for a particular case of a 2 Mq 
star. In order to illustrate the physics behind the existence of the 
loop at the end of the PMS we then discuss the limit case of a spe- 
cific value of the extent of the overshoot and the evolution of the 
chemical abundances that produce the loop. The paper ends by ad- 
dressing how different input physics and stellar parameters change 
the limit value of the overshoot producing the loop and what ob- 
servational evidence could be expected in order to identify a star 
undergoing such a phase in its PMS evolution. 



2 MODELLING CONVECTIVE OVERSHOOTING IN 
PMS EVOLUTION 

In this work we propose to study in detail the effects of convective 
overshooting in the PMS evolution of intermediate mass stars, as 
they approach the ZAMS, and so it is relevant to start by address- 
ing how overshooting is implemented in an evolutionary code. We 
will restrict our discussion to the case of intermediate mass stars 
(around 2 Mq), where there is a small central convective core. 

The boundary of a convective core Tco, with a mass M^o 
(see Fig. 0, is defined by the instability criteria, namely the 
Schwarzschild criterium or the Ledoux criterium if the chemi- 
cal composition is not homogeneous. For the Schwarzschild cri- 
terium the border between the convective zone (core) and the ra- 
diative zone (envelope) is located where the adiabatic gradient Vad 
equals the radiative gradient Viad- Convective elements that reach 
the boundary with non-zero velocities overshoot it, reaching a re- 
gion where Vrad<Vad. To transport all the energy generated in the 
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Figure 1. Schematic representation of the behaviour of temperature gradi- 
ent (thick continuous line) in the convective core and overshoot region of 
an intermediate mass star. The notation is defined in the text. CZ represents 
the proper convective zone, while MZ is the mixing zone. 

convective core upwards, the luminosity transported by radiation in 
the zone immediately above the boundary must be superior to the 
total luminosity, that is, in this zone one must have V rad<Vov<Vad . 

All mod els of overshooting have inadequacies teenzini" 19871 
ICanutdll997t) . as it is not possible to fully describe in one dimen- 
sion the effect of overshooting plumes. The standard option used is 
to parametrize the effect overshoot may have in the local structure 
by adopting a prescription. There are two major aspects of over- 
shoot that we need to consider in such a prescription if we want 
to discuss the effect on stellar evolution. These are the temperature 
stratification within the overshoot layer and the mixing taking place 
in such a region at the top of the "proper" convective core. 

2.1 Thermal stratification at the overshoot layer 

The temperature stratification below rco is given by the MLT for- 
mulation (which corresponds to the convective zone; CZ). Given 
the very small superadiabaticity expected in the proper convection 
zone the actual structure is not sensitive to the details of the par- 
ticular formulation we use, or to the mixing lengt h parameter a. 
Here we will adopt a value of a=1.35. As shown bv lMargues et alJ 
(2004) for stars with masses around 2 Mq any other reasonable 
value (as the one obtained from a solar calibration) will produce 
the same evolutionary tracks near the end of the PMS evolution. 

As discussed in several works, the temperature stratification 
in the overshoot layer i s expe cted to be subadiabatic. We assume, 
as in iDemargue et alj j20o3), that the stratification in the over- 
shoot layer is unaffecte d by the overshoot, a result also obtained by 
iBrowning et alJ i2004h . Consequently, within the overshoot layer 
of thickness dov (corresponding to rco^r^rco+dov) the stratifica- 
tion is assumed to be radiative (Vov=V,ad). Although it is a crude 
choice, its implementation in an evolutionary code is simple and 
does not affect the key aspects of overshoot that can be expected 
to interfere with the nuclear reactions that take place as the star 
approaches the main sequence. To confirm if this is so we also con- 
sider here other options for the thermal stratification of the over- 
shoot layer. 

The representation of overshooting should be implemented 
using more sophisti cated non-local formulations or the results of 
3D simulations (e.g. lGrossmanll995lCanutoll997tlBrowning et alJ 
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l2004b . Unfortunately, most of these formulations can not be fully 
implemented in detailed evolutionary codes as the non-local char- 
acter of overshooting compromises the efficiency of the integration 
in time for stellar evolution. 



Table 1. Mean lifetimes (in years) for the nuclei involved in the CNO bi- 
cycle, for X=0.7 and p=64 g cm~^ (from Icaughlan & Fowler 1963). 
Mean lifetimes for secondary nuclei are not included since their /3-decay 
is very fast (10^ - 10'^ s) compared to the mean lifetimes of the primary 
nuclei (T(i=T/10^ K). 



2.2 Mixing in the overshoot layer 

The key assumption, that we adopt here, is that the mixing of the 
stellar matter within the overshoot layer is highly efficient taking 
place at the same tim e scales as in the proper convection zone (see 
iBrowning e t al .'2004 V That implies that in terms of evolution, and 
in particular at the end of the PMS, it is fully justifiable to consider 
that convection and overshooting induce instantaneous mixing as 
far as evolution is concerned. Also, if the stellar mass is not too 
high, as is the case for the regime of intermediate mass stars we 
are addressing here, such mixing takes place in time scales much 
shorter than the mean lifetimes of all reactions in the CNO bi-cycle 
that may take place before the star reaches the ZAMS. 

As long as the assumption of instantaneous mixing within the 
overshoot layer stands, the implications for the evolution of more 
complex formulations for the thermal stratification of overshoot are 
equivalent to the results provided by the simpler formulation that 
is being adopted. For the phase of evolution we discuss here the 
classic assumption of using the limit of instantaneous mixing is ad- 
equate as it provides a representation of a lower limit for the extent 
of overshoot that may produce the loop. This extent is well below 
the ex pected extent of overshoot for this mass range iRibas et'all 
I2OOOI) . 

In order to establish to what extent our results depended on 
the mixing model in the overshoot r egion, we also consider models 
calculated using the prescription of IVentura et alJ il998h . Accord- 
ingly, the evolution of element i follows the diffusion equation: 
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The diffusion coefficient D is approximated in a convective 
zone by D=u Zd/3, where u is the average turbulent velocit y, com - 
puted according to eqs (88), (89) and (90) of Canuto et al1 <199d) . 
and /(J is the convective scale length. The convective scale length U 
is given by 



(2) 



where Zup is the distance from the top of the convective zone in- 



creased by (3H^p and analogously for z\o 



. Here, /3 is a fine tuning 



parameter, necessary for an exact fit of actual stars, and Hp the 
pressure scale height at the top of the convective region. The pa- 
rameter P is constrained by /3<0.25. 

For diffusive overshooting, we write: 



u — Uh exp ± 



C /thick 



In 



f- 



(3) 



where Uh and Pb are the turbulent velocity and the pressure at the 
border of the convective zone, P is the local pressure, is a free 
parameter and /thick is the thickness of the convective region in frac- 
tions of the local Hp. This way, the turbulent velocity decreases 
exponentially outside convective regions. The parameter C_ controls 
the e-folding distance with which the velocity of convective eddies 
decay outside convective regions; a higher C, means that the veloc- 
ity of the convective eddies decays slower and therefore a bigger 
region is affected by partial convective mixing. The diffusive scale 
is approximated by li=j3Hp in overshoot regions. 
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2.3 Extent of the overshoot layer 

Following the assumptions discussed above the only aspect remain- 
ing to be defined is the extent dov of the overshoot l ayer. We adopt 
the co mm on prescription (s ee Maedei 1975; M owlavi & Forestinil 
I1994I and iDemargue et ant2004J that parametrises the extent of 
the overshooting layer in units of the local pressure scale height 
Hp=\Ar /d In P| . So, the extension of the overshooting is given by 
dov=Q!ov-Min {Hp, no), where the border r^o of the "proper" con- 
vective core corresponds to the radius defined by the Schwarzschild 
criterium and ciov is a free parameter to be defined. We note, how- 
ever, that eclipsing binarie s data indicates t hat this parameter is 
dependent on stellar mass (Ribas et al. "200^, reaching a value of 
about 0.2 for A/~2 A/q. As discussed above, inside rniz='"co+dov 
there is complete mixing (we shall call this zone the mixed zone - 
MZ). 



3 OVERSHOOTING AND EVOLUTION OF THE 
INTERNAL STRUCTURE 

Models of a 2 A/0 star with values of varying from to 
0.13 have been produced (see Table|2j. The evolution was calcu- 
lated using an initial helium abundance yo=0.28, initial metallicity 
2'o=0.02 and mixing length parameter a=1.35. Mo dels were cal - 
culated using the CESAM stellar evolutionary code jMoreil99"ii . 
The up-to-date physical ingredients and adopted numerical proce- 
dures relevant for the early evolution near the ZAMS - used to build 
the models discussed here - are described in detail bv lMargues et alJ 
C2004) . 

The key aspect in the calculation of models near the ZAMS 
is the energy production. For the mass range discussed here this is 
the CNO bi-cycle which consists of two cycles. In the first, the CN 
cycle, only carbon and nitrogen are involved: 
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The second, the ON cycle, is entered by the termination of the last 
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Figure 2. Left panels: evolution of the convective zones for models with different values of aov'. aov=0 (top panels), Oov=0.1295 (middle panels) and 
aov=0.1300 (bottom panels). Cross-hatched areas represent the extent of overshooting. Right panels: evolutionary tracks in the HR diagram for the same 
cases as in the left panels. Symbols on the track correspond to evolutionary ages of special interest: they ai'e also shown in the left panels (see Tablel2l. 



reaction of ^4) through the 7 channel (instead of the a channel): 

^^N + 'H ^ "0 + 7 
'^O + 'H ^ '''F + 7 

^^O + ^H ^ "N + *He. (5) 

TableQshows, however, that the termination of the fusion of ^^N 
through the a channel is much more likely (about 10* times) than 
through the 7 channel. The equilibrium timescale of the CN cycle is 
therefore much shorter than that of the full CNO cycle; for the PMS 
phase, only the processes in ^4) - the CN cycle - will be considered 
in our discussion (although the models use the full network). There 
is still another ON cycle, which is entered by the termination of 
the last reaction of ^5} (the fusion of ^^O) through the 7 channel; 
its importance, however, is very small. TableQalso shows that the 
slowest process in ^4) is, by far, the combustion of ^*N. Therefore, 
before the abundances in CN elements reach equilibrium almost all 
original ^"^C and ^'^C must be burned into ^*N. 

All models discussed here use the full network of reactions 
(PP chains and CNO cycle) to calculate the energy production. The 
PMS birthline (Palla & Stabler 1991) is not used in these calcula- 
tions as the chemical profile at the core of the models for a 2 Mq 



Table 2. Age (in Myrs) of the models at different phases of evolution for 
different values of Qov. These phases correspond to: a - the central MZ 
appears; b - the central MZ reaches the first maximum; c - minimum extent 
of the central MZ; d - ZAMS. Also indicated is the type of symbol used in 
Fig.|2|to locate these phases. 



Phase Cfov Symbol 

0.0000 0.1295 0.1300 (Fig.|2) 
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9.2 


black circle 






9.7 
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black diamond 
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12.1 


star 



star near the ZAMS are not significantly changed by what happens 
at the very early stages of stellar evolution. 

Figure |2| shows evolutionary tracks for several values of Oov 
The most remarkable feature is that there are two kinds of tracks, 
with or without a "loop" just before the ZAMS. The transition be- 
tween these two types of tracks happens at a very definite value of 
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Figure 3. Evolution of the production of energy (top panel) and the central 
abundances (lower panel) without overshooting. 



Oov: for q:ov=0.1295 there is no "loop" in the evolutionary tracks 
on the HR diagram, while there is a "loop" for aov=0.1300. For 
values of Qov lower than 0.1295, the tracks do not differ signifi- 
cantly; the same happens for values of Oov higher than 0. 1300. This 
lo op can also be seen in the evolutionary models available online, 
bv lSiess etafj<200a) , for a 2 Mq star with Z=0.02 and a„y=0.2. 
This feature requires an explanation. In the following sections we 
will describe in detail the evolution of the internal structure of a 
2 Mq star without and with overshooting in order to understand 
the origin of these features in the models of PMS evolution. 



3.1 Evolution without overshooting 

The evolutionary track of a 2 Mq star on the HR diagram without 
overshooting is shown in Fig.|2|(upper panel), as well as the evolu- 
tion of the convective zones (top right panel), during the last phase 
of the PMS. 

At an age t^,=5 Myrs, the star is completely radiative and all 
its layers are contracting. The central temperature rises, until it be- 
comes high enough to bum ^He and ^^C. The combustion of both 
^He and ^^C depends strongly on the temperature (e^T^*^ to T^**); 
therefore, the energy generated is strongly concentrated in the cen- 
tre of the star. The energy flux from the central regions becomes 
so high that it can not be transported by radiation alone. A central 
convective zone appears at this point. Figure|3]shows the evolution 
of the production of energy and the central abundances of CN ele- 
ments. Comparing Figs.|2|(top left panel) and|3| it is clear that the 
convective core appears at the time when ^^C-buming starts (at an 
age ti,—6.6 Myrs); this instant is indicated in Fig.|2|by the black 
triangle (see Table|2}. 

The structure of the core implies that Pc/p is about 6 for a 
convective model, while being higher than 50 for a radiative one. 
Because the convective model has to be less centrally concentrated, 
the central zone must expand when convection appears; this expan- 
sion absorbs practically all the energy produced in the centre by 



^^C-buming, which does not, thus, reach the surface. We have the 
apparently paradoxical result that when nuclear energy sources be- 
come important the total luminosity of the star decreases (see top 
right panel of Fig.|2j, although there is now a new energy source. 

As the burning of ^^C becomes more efficient with rising cen- 
tral temperature and density, the central convective zone grows; 
when the abundance of ^^C in this zone decreases, the energy gen- 
erated in the centre of the star decreases as well. In Fig.|2| the black 
square indicates the instant the central convective zone reaches its 
maximum extent (t*=7.2 Myrs). As the energy generated in the 
centre of the star decreases, the extent of the central convective 
zone decreases. As discussed previously, when this happens the re- 
gion of the convective zone must now strongly contract in order to 
adapt to the structure of a radiative core. The central density rises 
again, as well as the central temperature. 

For a 1 A'Iq star, PP chains would generate now practically 
all the energy, since ^''N-burning would be always too slow for the 
CN cycle J4j to generate a large fraction of the energy (see TableQ 
for T(j — 15). For a 2 Mq star, however, PP chains can not stop the 
contraction of the stellar core. As almost all ^^C originally in the 
core has been burned and the convective core retracts as a conse- 
quence, the central zone contracts again. The luminosity generated 
by this contraction causes the luminosity of the star to increase, as 
can be seen in the top right panel of Fig. |2|and in Fig.|3| As the 
central temperature rises, the combustion of becomes more 
efficient (Table 0; again, the energy generation becomes concen- 
trated in the centre of the star; the convective zone grows and, as a 
consequence, the luminosity of the star decreases again. The instant 
when the convective zone reaches a minimum is indicated in Fig.|2| 
by the black circle (t*=8.7 Myrs). The contribution of the CNO 
cycle J4j to the energy generation grows, until the star is finally in 
hydrostatic equilibrium. The ZAMS is indicated in Fig.|2| top right 
panel, by the star (t* = 10.1 Myrs). 

Figure |3| shows that the abundances of ^^C and ^"'C do not 
grow significantly while the ^*N-burning reaction is becoming 
more efficient; all ^^C produced by the burning of ^''N very rapidly 
burns into ^^C, which even more rapidly burns into ^''N again (Ta- 
ble0. The burning of ^''N, because it is so slow compared with the 
other reactions of the CN cycle J4j, acts as a bottleneck. 



3.2 Evolution with overshooting 

As seen above, for aov^O.1300 there is a loop in the evolu- 
tionary track in the HR diagram during the final stages of the 
PMS. For Qov ^0-1295 there is no such loop; the transition is 
quite abrupt. Figure|2|shows the evolution of the convective zones 
with Qov=0.1295 (middle panels); it is qualitatively similar to the 
evolution with aoy—0. The two main differences are that with 
Qov=0.1295 the extent of the MZ is higher and the evolutionary 
time is longer (see Table|2j- The extent of the MZ at the minimum 
after the main ^^C-buming is also significantly higher. These dif- 
ferences are easily explained by the extra extent of the MZ caused 
by the overshooting; the burning of ^^C takes more time because 
the MZ extends further and therefore there is more ^^C to burn, 
while the rate of burning is similar with aov=0-0 and aov=0.1295 
(since the central densities and temperatures are similar). Table |2| 
shows the ages of the models at some stages of the evolution. Be- 
fore the appearance of the central MZ, the evolutionary times are 
the same. 

Figure|2|shows the evolution in the HR diagram (bottom right 
panel) and the evolution of the convective zones (bottom left panel) 
of a model with aov=0.1300. In the bottom right panel of Fig.|2| 
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Figure 4. Evolution of the production of energy (top panel) and central 
abundances (lower panel) with aov=0.1300. 



the black triangle shows the instant the central MZ appears for 
the first time {ti,=6.& Myrs); this zone reaches its maximum ex- 
tent around t*=7.2 Myrs (black square). As the ^^C present in 
the central MZ bums out, the MZ retreats, reaching a minimum 
at i*=9.2 Myrs (black circle). It grows again as the burning of 
becomes more efficient, and, while with aov=0.1295 the reactions 
of the CNO cycle J4} reach equilibrium and the central MZ sta- 
bilises, with aov=0.1300 the MZ grows suddenly. Until this point, 
the evolution for Qov=0.1295 and aov=0.1300 is the same; here, 
they split completely. This event (at t*=9.7 Myrs) is indicated in 
Fig. |2| with the open square. The central MZ grows rapidly, reach- 
ing a maximum (~0.7 M0=O.35 !) at f*=9.8 Myrs (black di- 
amond). This sudden growth of the MZ is caused by an increase in 
the central abundance of ^^C and ^^C, which are burned very effi- 
ciently because the central temperature is now higher than during 
the main ^^C-burning phase; so, the energy produced in the centre 
of the star increases rapidly. 

Figure|4|shows the evolution of the production of energy and 
the central abundances with aov=0.1300; the main difference to 
Fig.|3|(cfov=0) is the sudden rise in the abundance in ^^C and ^'^C 
(as well as a drop in the abundance in ^*N). This causes the rise 
in the production of energy through the CNO cycle seen in the left 
panel of Fig.|4| causing in turn the growth of the central MZ seen in 
Fig.m From this point, the MZ retracts again, reaching a new mini- 
mum at f* = 10.5 Myrs (black pentagon). The model of the star then 
returns to the "normal" track, reaching the ZAMS at =12.1 Myrs. 

The MZ retracts so fast because it grew so fast before. This 
sudden growth causes an expansion of the central regions, which 
leads to a big drop in the central temperature and particularly a 
big drop in the central density. The reaction rates of ^^C and ^'^C- 
burning, depending so strongly on the temperature and density, 
drop rapidly. Since this sequence of events (after the sudden growth 
of the central abundance in ^^C and ^"^C) does not depend on the 
extent of the central MZ when aov^O.1300, depending only on the 
central temperatures and densities, the format of the loop on the 



HR diagram does not depend on the actual extent of overshooting 
above aov=0.1300. 

3.3 The central abundances 

In this section we study the evolution of the central abundances 
for Qov=0 and Qov=0.1295. The value aov=0.1295 is particularly 
interesting because it is almost enough to produce a "loop" in the 
evolutionary tracks on the HR diagram. The comparison between 
the cases aov=0 and q:ov=0.1295 shows the reasons why a "loop" 
is produced for q:ovSs0.1300 and why is this transition so abrupt. 

The chemical composition of ^^C at the central convec- 
tive zone reaches a minimum at f*=8.7 Myrs for otm=G and at 
U=Q.l Myrs for aov=0.1295. The ^^C is not fully burned above 
the MZ at this stage. If temperatures above the MZ remained too 
low to burn significant amounts of ^^C there during the final expan- 
sion of the MZ, the MZ would add this ^^C during the following 
expansion. Figure|5|(upper panels) shows the evolution of the tem- 
perature of the BMZ (labelled Tbmz) for Qov=0, 0.1295, 0.13. The 
temperature (Tign) at which the mean lifetime for ^^C is 0.75 Myrs 
(for the density at the BMZ), is also shown. It is clear that the tem- 
perature at the BMZ decreases as the amount of overshooting in- 
creases; overshooting extends the MZ, so the BMZ is farther from 
the centre of the star, where the temperature is lower. Without over- 
shooting, the temperature of the BMZ is higher (during the final 
expansion of the MZ) than the temperature needed to significantly 
deplete ^^C before the MZ adds it. So we expect that the MZ adds 
almost no ^"^C during its final expansion. 

As we increase the amount of overshooting, the temperature 
at the BMZ decreases (see upper panels of Fig.|5), making ^^C- 
buming above the MZ incomplete. As a result, some ^^C is added 
to the MZ during its expansion. Finally, when q:ov=0.1295 there is 
almost no depletion of ^^C above the MZ during its last expansion. 
The increase with Qov of the amount of ^^C added to the MZ when 
Qov approaches the critical value of 0.I295 is very high. Even a 
small increase of Qov from Qov=0.1290 to Qfov=0.1295 is enough 
to increase noticeably the amount of ^^C added to the central MZ. 

If there is no ^^C added to the core (as it is the case with- 
out overshooting), the CNO reactions reach equilibrium shortly af- 
ter ^*N-burning becomes efficient enough (that is, shortly after the 
central MZ starts to grow for the last time). When, due to the extra 
extent of the central MZ caused by overshooting there is more ^^C 
added to the MZ than that produced by ^^N-buming, there must be 
an excess of ^'^C (and '^■'C)-burning to keep the relative abundances 
of CN elements at equilibrium. 

This excess burning of ^^C (and ^'^C, since the extra burn- 
ing of ^^C generates an excess in the abundance of ^''C relative to 
the equilibrium abundances) causes an increase on the luminosity 
produced within the central MZ; to transport this extra luminosity 
out, the central MZ must grow. Figure|S|shows that only when Qov 
approaches the critical value does this increase in the luminosity 
becomes noticeable. The growth of the central MZ must bring in 
turn an even greater amount of ^^C to the MZ, causing an even 
greater increase in the luminosity produced within the central MZ. 
This is only stopped because such a fast increase in the energy pro- 
duced within the core causes it to expand and cool, making CNO 
reactions much less efficient. 

In short, the abruptness of the transition from an evolution 
with a loop to an evolution without a loop comes from the sensibil- 
ity of the CNO reactions to the temperature, as the parameter Qov 
determines the temperature of the BMZ, which in turn determines 
how much ^^C is burned before the MZ grows. 
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Figure 5. Upper panels: evolution of the temperature of the BMZ (full line) for several values of aov (where T^=T/1Q^). The temperature (Tign) at which 
the mean Ufetime for ^-^C is 0.75 Myrs (for the density at the BMZ), is also shown as a dashed line. Lower panels: the abundance of ^^C the BMZ finds 
during the final expansion of the MZ (dashed line); evolution of the central abundance in ^^C during the same period (full line); and the abundance in ^^C 
the BMZ would find if there would be no depletion of '^^C above the MZ due to nuclear reactions (dotted line). The shadow in the upper panels indicates the 
time range of the corresponding lower panel. 



4 SENSITIVITY OF THE RESULTS TO THE INPUT 
PHYSICS 

The physics of the early evolution for a 2 Mq star is expected to 
be simpler than for the more advanced stages of evolution. Aspects 
as diffusion and settling as well as the existence of steep chemical 
gradients as a result of hydrogen burning in the main sequence are 
not important in the cases discussed here. In this section we address 
some of the key aspects that may be relevant and the possible effect 
these may have on the existence of the loop in real stars. 

Although we presented a very definite value of ciov beyond 
which a loop is produced in the evolutionary tracks on the HR di- 
agram, this "critical" value can not be taken at face value, since 
it depends on several known and unknown factors. Such a critical 
value depends mainly on the nuclear reaction rates, but other as- 
pects of the models may modify the actual number. We stress that 
the main goal of this work is not the determination of a critical 
value for Oov but to establish what are the physical and numeri- 
cal ingredients responsible for producing the loop in the evolution 
calculation. 



4.1 Physics and numerics of the overshoot layer 

Numerical simulations and observational tests (see Section^ con- 
sistently indicate that stars are expected to have mixing zones that 
go beyond the border predicted by an instability criterium for con- 
vection. The unresolved questions are the extent and the physics 
of these extra mixing zones. Here we have included such an ex- 
tra mixing zone (overshoot layer) at the top of a small convective 
core by adding an extra mixing layer of size aov Hp to the proper 
convection zone modelled according to the mixing length theory. 

There are different prescriptions for overshooting as discussed 
in Section |2| which will provide slightly different results for the 
same value of a^w In order to confirm that our results are not af- 
fected by the thermal stratification adopted within the overshoot 
layer we have also calculated the evolution when the limit of 
Vov=Vad is used within the overshoot region. As long as instanta- 
neous mixing is assumed the critical value of ciov does not change 
by more than about 0.0005. The implications for the evolutionary 
track in the HR diagram are not changed and the existence of the 
loop is still confirmed in this limit. The actual stratification in the 
overshoot layer of stars is expected to lie between these two val- 
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Figure 6. Evolution of the luminosity produced within the central MZ by 
the burning of (full line), (dashed line) and (dotted line). The 
panels coiTespond to different values of Qov ■ 

ues (Vrad<Vov< Vad). Consequently, any temperature profile will 
produce tiie same beiiaviour we discuss in tiiis work. 

Models calculated witii different precision (number of sliells 
for ea ch model and time step of the evolution - see Mar ques et alj 
12004 ) or using different numerical methods yield different values 
for the critical Ofov These differences (below 0.001) are a direct 
consequence of the need for a careful numerical scheme to follow 
in time and space the border of convective and mixing zones in stel- 
lar interiors. These numerical tests indicate that the existence of the 
loop is to be expected when adequate numeric procedures are used 
for time and space integration. As far as this work is concerned 
these uncertainties are small corrections which will only become 
relevant if the measured extent of overshoot in stars is close to the 
critical value. As discussed above (Section 0, observational con- 
straints on the value of the extent of an overshoot layer are not 
clear yet. Some works dRibas et all200d) indicate that the amount 
of overshooting in the cores of 2 Mq stars can be significantly 
larger than aov~0.13. 

4.2 Diffusive versus instantaneous mixing 

Here we have adopted instantaneous mixing for the overshoot layer. 
This is a crucial assumption that needs to be confirmed. 

The time scale of the relevant nuclear reactions is signifi- 
cantly larger than the expected convective turnover time for the 
core. However, the mixing at the overshoot layer, in particular at 
the very top, may have longer mixing timescales when compared 
with the proper convec tive zone. Models in corporating overshoot 
as a diffusion process <VenturaetJ]ll998h indicate that for stars 
with convective cores the effect on the mixing is equivalent to an 
instantaneous overshooting around 0.2 Hp. This result has been 



reinforced bv I Ventura & Castellanil i2005h where nuclear burning 
and mixing were self-consistently coupled for post-main sequence 
stars and compared with instantaneous mixing. The confirmation 
that such a behaviour is also valid for the pre-main sequence evo- 
lution is required, but given the physics at play in this early phase 
no major differences are expected. 

In order to verify this assumption we have also considered the 
case when mixing is not instantaneous outside convective zones. 
However, the exponential decay of the velocities outside the con- 
vective zone means that over time the mixed zone is larger. During 
the first retreat of the central CZ extra ^^C is brought to the core 
so the energy generation there is higher than it would if no over- 
shooting were present. That keeps the convective zone bigger and 
the temperature at its border lower, so when the CZ grows again it 
finds above a significant abundance of ^^C, as happens in the case 
of instantaneous mixing in the overshoot zone. 

Consequently, a loop is also produce d on the evolu t ionary 
tracks calculated using the prescription of Iventura et alJ il998h 
with values of (^JsO.Ol and /3^0.05. These values shou ld be com- 
pared to C=0.03 used bv IVentura & Castellanil <2005h to fit the 
cluster NGC 1866. The loop is found to be very similar to the one 
obtained with instantaneous mixing; both the trajectory in the HR 
diagram and the time it takes. 

4.3 Stellar rotation 

Rotation can also play a role on how overshoot behaves. Consider- 
ing the aim of this work, the important questions to be addressed are 
the following: can rotation change, considerably, the time scale of 
mixing when compared to the burning time scale? Is the efficiency 
of m ixing affected by the rotation? 

iBrowning e t al.' 1 2004) presented three-dimensional simula- 
tions of core convection for a 2 Mq star. They found that the core 
has a differential rotation and an overturning period of the global 
scale convection of about one month. On the other hand they sup- 
port that overshooting and penetrative convection are both effec- 
tive in mixing the chemical composition. The existence of an ef- 
fic ient overshootin g layer for rotating stellar cores was also found 
bv lDeuDredll998h using 2D simulations for a 8.75 Mq star. Con- 
sequently there is no evidence to support (on the contrary) that ro- 
tation will affect the results reported here, and in particularly our 
assumption concerning the in stantaneous mixing . 

It is worth noticing that lMevnet & MaedeJ i200(jl) found for 
high mass stars (9^M/M0^12O) that rotation modifies the tracks 
in the HR diagram as moderate overshoot would do. A similar con- 
clusion is reached by Noels et al. (2004) for A stars, where the 
authors argue that rotation reduces the efficiency of convection 
and consequently the extent of overshoo ting. However the obser- 
vational evidence jHerbst & Mundl2005h indicates that most PMS 
stars loose their angular momentum well before arriving at the main 
sequence, when they are already slow rotators. In such cases sup- 
pression of overshooting due to rotation is weaker and so a loop 
may be expected in intermediate mass PMS slow rotators close to 
the main sequence. 



4.4 Dependence on stellar parameters 

The critical value of the parameter Qov depends strongly on the stel- 
lar mass (or the size of the convective core), decreasing for higher 
stellar masses. For a 3 Mq star, there is a loop in the evolution- 
ary track on the HR diagram for Qov ^0.0593, while for a 4 Mq 
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star the same thing happens for aov^O.0391. So, for higher mass 
stars a relatively small amount of overshooting is enough to pro- 
duce a loop before the ZAMS. This is due to the faster evolution of 
higher mass stars, leaving less time for the depletion of ^^C in the 
regions above the MZ before the MZ grows again. So Tign, which 
we defined as the temperature at which the mean lifetime of ^^C IS 
Tci2=0.75 Myrs, should be replaced by the temperature at which 
Tci2=0.27 Myrs for a 3 Mq star and rci2=0.13 Myrs for a 4 Mq 
star. The region where the temperature is equal to Tign lies there- 
fore nearer the centre of the star, and a relatively small amount of 
overshooting is enough to bring the BMZ to those regions. 

The evolution is faster for higher masses; the loop lasts about 
4x10^ yrs for a 3 Mq star and about 2x10^ yrs for a 4 Mq 
star, making a possible observation much more difficult. An ad- 
diti onal difficulty for the case of a 4 Mq star is that the birthline 
(see lPalla & StahleJll99lLll992h is closer to the ZAMS; the loop 
should happen at an age of about 5x10^ yrs after the birth of the 
star. The effects of the protostar phase could be of some impor- 
tance at such an early age. Also, the assumption of instantaneous 
mixing becomes weaker as we move towards higher mass models. 
In such cases the requirement to consistently follow the mixing and 
the nuclear reactions becomes more important, being necessary to 
calculate models wi th time dependent mixing co upled with nuclear 
burning (as done in I Ventura & CastellanilEoOS . for post-main se- 
quence stars) at the central core of these young stars. 



5 EXPECTED OBSERVATIONAL EVIDENCES 

The observational confirmation that a particular star is in the loop 
can be straight forward, but the probability of finding such a star 
is very small as it depends on the relative time that a star spends 
in the loop. Very young clusters with stars in the range of 1.5 to 
3.0 Mq arriving at the ZAMS are preferential targets for finding 
loop stars. The location at the HR diagram (see Fig. has to be 
complemented with the determination of the mass of the star. Con- 
sequently the study of possible candidates in binary systems, with 
a smaller mass component to determine the age precisely, is the 
natural observational test for the existence and determination of the 
characteristics of loop-stars. 

An alternative to study binaries, is to obtain seismic data on 
several individual stars in order to discriminate between stars with 
the same radius but different masses. Such a test can be applied 
as long as a few oscillations can be identified in order to provide 
the observational value for the large frequency separation At for 
modes of degree £ and calculated at a referenc e frequency (e.g. 
iMonteirollioollFemandes & Monteirdl2003^ . This quantity (for 
^=0) can differ by as much as 6 fiHz between a loop-star and a 
star without overshoot at the same location in the HR diagram. 

If very precise seismic data are available the small frequency 
separation, can also be determined using modes of degree 

£ and £+2 (being calculated at the same reference frequency). The 
small freque ncy separation, being sensitive to the central regions of 
the star (e.g. lMonteiro et aOboOa) . can differ by as much as 30% 
at i/r=1200 fiHz for location 2 in Fig.0 providing a powerful ad- 
ditional test to the large frequency separation. 

Another seismic indicator, mainly sensitive to the central re- 
gions of the stars (Roxburgh & Vorontsov 2003; Roxburgh 2005), 
is the ratio ro2=5o2/Ai of the small separation (using modes of 
degree and 2) and the large separation (for modes of degree 1). 
This quantity can provide an even more definite discriminant for 
stars with strong structural differences in the interior but with the 




Figure 7. Evolutionary tracks in the HR diagram for thi'ee stars. One with 
M=2 Mq (full line) having a loop due to the presence of overshoot, and 
two stars of smaller mass (dashed lines) without overshoot which cross the 
track of the more massive star. The point "0" marks the beginning of the 
loop, the point "4" (in the same location on the HR diagram) marks the end. 



same radius. This is the case between loop-stars and the correspon- 
dent star without overshoot at the same location in the HR diagram. 
The values of ro2 for models occupying the same location in the HR 
diagram show differences that can be as high as 0.025 (~26%). 



6 CONCLUSIONS 

We have shown that a moderate amount of overshooting in evolu- 
tion models of young intermediate mass stars can cause a loop in 
the evolutionary tracks on the HR diagram just before the ZAMS. 
This is a consequence of using a more precise numerical scheme for 
calculating the time-step together with the best up-to-date physics. 
The loop corresponds, for a 2 Mq star, to a drop of about 35% 
in luminosity and about 1000 K in effective temperature; the loop 
lasts about 1.5 Myrs. A detailed description on why the full nuclear 
network predicts in the models such a loop has been presented. 
This analysis indicates that such a behaviour may be expected if 
the standard formulations for core overshoot in stellar models are 
representative of the average role convective overshooting has in 
the early evolution of intermediate mass stars. The major effect of 
the existence of the loop is on the age of the star at the ZAMS, be- 
ing higher by as much as 20% relative to the evolution without the 
loop. 

The observational identification of loop-stars can provide def- 
inite tests on some of the key aspects of the physics that determine 
the evolution in early and main sequence intermediate mass stars. In 
particular the existence and extent of overshoot at the core of these 
stars and the reaction rates for the branches in the CNO cycle re- 
sponsible for producing the loop. Such a combination of precise ob- 
servational constraints make the existence of loop-stars an interest- 
ing and worth pursuing observational opportunity to test the mod- 
elling of stellar structure and evolution at these early stages. The 
observational challenge is the short duration of the loop; less than 
about 2 Myrs. The loop is even shorter for higher stellar masses, 
for which the amount of overshooting needed to produce a loop be- 
comes smaller and therefore more likely. Further analysis on how 
overshooting and nuclear reactions are coupled in the higher mass 
models is required in order to identify if the loop can be expected 
in this regime. 
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With the forth coming a sterose ismic space missions (e.g. 
iBaglin & and the Co RoT Teami |2003ll and ongoing ground based 
campaigns for the seismic study of stars across the HR diagram, 
there will be a large set of young stars whose seismic properties 
will be determined with sufficiently high precision to be possible 
to use the seismic analysis we have discussed here. There are in 
particular campaigns planned for the seismic observation of young 
clusters whose ages are estimated to be around the expected values 
for the existence of loop stars (from about 5 to 12 Myrs). Con- 
sequently, the possibility of confirming the existence of loop stars 
may be possible in the not so distant future. The implications of 
finding a few such stars in a young cluster will be of great rele- 
vance for the modelling of the early evolution of intermediate mass 
stars. 
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